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需求 不 确定 的 电动 汽车 换 电站 选 址 鲁 棒 模型 ， 
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摘 要 : 在 电动 汽车 换 电站 选 址 问题 中 ， 用 户 充电 需求 是 重要 的 输入 参数 ， 而 需求 受到 技术 发 展 、 电 动 汽车 保有 量 、 
国家 相关 政策 等 因素 的 影响 ， 通 常 较 难 准 确 预 测 。 本 文 假设 交通 网 络 上 用 户 需 求 是 不 确定 的 ， 利 用 基 约 束 鲁 棒 优化 方 
法 提出 基于 用 户 路 径流 量 不 确定 的 换 电 站 选 址 鲁 棒 模型 ， 并 将 鲁 棒 模型 转化 为 等 价 的 线性 规划 。 最 后 应 用 
Nguyen-Dupius 网 络 作 为 数值 算 例 进行 分 析 ， 结 果 表 明 当 解 的 鲁 棒 水 平 较 低 时 ， 适 当 增 加 选 址 成 本 ， 可 以 较 大 程度 
提高 解 的 鲁 棒 性 。 同 时 ， uou de e 随 着 鲁 棒 水 平 的 提高 ， 换 电站 选 址 成 本 对 路 径流 量 的 波动 
越 来 越 敏感 ， 续 航 里 程 的 改变 对 换 电 站 选 址 策略 影响 较 大 。 
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Research on robust battery swap stations location model of electric vehicles 
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Abstract: The users' charging demand is important input parameters in the battery swap stations location problem of electric 


vehicles. The demand is influenced by technological development, the number of electric vehicles and the relevant policies, 


and it is often difficult to predict accurately. This research examines the robust battery swap stations location problem with 


es flow uncertainty. Supposing the users' demand on the traffic network path is uncertain, the cardinality constrained robust 
approach is employed to formulate the robust battery swap stations location model. The model can be transformed to a linear 
optimization problem. We use the Nguyen-Dupius network as a numerical example. The results show that when the robust 
level of the solution is low, the robust level can be improved greatly by appropriately increasing the location cost. In addition, 
the sensitivity analysis is conducted for important parameters of the model. The results show that, the location cost is 
increasingly sensitive to the fluctuation of path flow with the improvement of robust level. The change of driving range has a 
great impact on the location strategy. 
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电动 汽车 相对 传统 汽车 来 说 还 是 新 生 事 物 ， 在 实际 使 用 环节 中 
的 充电 设施 配套 体系 还 不 完善 ， 充 电 难 已 经 成 为 影响 消费 者 购 
有 E 动 汽车 因 其 清洁 节能 的 显著 优势 ， 逐 步 成 为 汽车 市 场 的 。 买 电动 汽车 的 最 大 障碍 帆 。 而 已 建成 的 部 分 充电 设施 却 呈 现 出 
主流 发 展 方向 。 在 巨大 的 环境 压力 和 未 来 市 场 空间 预期 的 双重 使 用 效率 低 甚至 闲置 的 状况 , 形成 了 当前 电动 汽车 发 展 的 “充电 
作用 下 ， 各 国政 府 纷 纷 介入 电动 汽车 市 场 ， 并 将 发 展 电动 汽车 。” 怪圈 ”外 ， 即 电动 汽车 充电 难 与 充电 设施 闲置 现象 并 存 。 因 此 ， 
作为 应 对 能 源 和 环境 问题 的 重要 举措 。 然 而 ， 电 动 汽车 市 场 的 充电 设施 的 选 址 布局 成 为 当前 吸 待 解决 的 问题 。 
推广 进展 却 很 缓慢 ,其 中 一 个 重要 原因 是 电动 汽车 续航 里 程 短 ， 换 电站 由 于 具有 种 种 优势 而 备 受 企业 和 用 户 的 喜爱 ， 如 换 
使 其 一 出 生 就 遭遇 了 里 程 焦虑 症 ( 指 驾驶 电动 汽车 时 因 担 心 突 电站 可 用 来 做 储 能 电站 、 削 峰 填 谷 、 电 网 谐 波 、 提 高 电网 负荷 
然 没 电 引起 的 精神 痛苦 和 忧虑 )。 传 统 燃油 汽车 也 有 续航 里 程 问 ” 率 、 可 对 电池 进行 慢 充 延长 电池 寿命 、 可 用 谷底 电价 进行 充电 
题 ， 只 不 过 目前 加 油 站 星罗棋布 ， 已 使 这 个 问题 得 以 解决 。 而 “降低 成 本 等 。 目 前 ， 国 内 建 换 电站 的 企业 主要 是 两 家 : 一 是 北 
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京 的 北汽 ， 二 是 浙江 的 时 空 。 北 汽 的 换 电 模 式 是 底盘 换 电 ， 主 
要 应 用 在 出 租车 上 ; 时 空 是 侧 方 换 电 , 主要 支撑 的 也 是 出 租车 、 
物流 车 等 工具 用 车 。 相 关 研 究 表明 在 电动 汽车 的 体 量 和 稳定 且 
迫切 的 能 量 补给 需求 下 ， 强 大 的 服务 能 力 很 快 能 让 换 电站 显现 
出 比 充电 站 更 强 的 熏 利 能 力 B。 然 而 ， 换 电站 投入 成 本 高 ， 建 
设 难 度 大 ,为 了 避免 电动 汽车 充电 难 与 充电 设施 闲置 并 存 的 “ 充 
怪圈 ”研究 如 何在 有 限 的 资金 和 紧张 的 土地 资源 下 , 合理 规 
划 布 局 电动 汽车 换 电站 至 关 重 要 。 

已 有 学 者 对 充电 设施 的 选 址 布局 问题 做 了 相关 研究 。 
Wang 只 假设 电动 旅行 车 的 行驶 路 径 已 知 ， 研 究 电池 换 电站 选 址 
优化 问题 。MirHassani 等 人 外 提出 更 为 灵活 的 线性 规划 模型 来 
解决 能 源 补充 设施 选 址 问题 。Jung 等 人 中 针对 电动 出 租车 充电 
站 的 布局 问题 提出 双 层 仿真 优化 模型 ， 该 模型 考虑 行驶 路 线 截 
流 问 题 和 排队 延迟 问题 。Dong 等 人 中研 究 基于 出 行者 活动 的 充 
设施 优化 布局 问题 ， 并 设计 遗传 算法 求解 模型 。Zheng 等 
[研究 续航 里 程 约束 下 的 电动 汽车 交通 流量 均衡 和 最 优 充电 设 
施 位 置 ， 提 出 双 层 模型 ， 上 层 最 优化 充电 设施 位 置 ， 下 层 是 流 
量 均衡 。Huang 等 人 号 考虑 交通 网 络 上 起 该 点 之 间 存 在 多 条 路 
径 ， 提 出 基于 流量 的 集 覆 盖 模型 。He 等 人 (9 对 北京 充电 设施 
的 选 址 布局 问题 做 案例 研究 ， 其 目标 是 将 当地 的 供应 限制 和 对 
充电 设施 的 需求 纳入 三 种 选 址 模型 ， 即 集 履 盖 问 题 、 最 大 履 盖 
问题 和 了 -中 位 问题 。 杨 表 等 人 5 提出 了 电动 汽车 物流 配送 系统 
的 换 电站 选 址 与 配送 路 径 优化 问题 ， 建 立 了 整数 规划 模型 。 陆 
坚毅 等 人 0 从 充电 站 运营 商 的 角度 出 发 ， 在 实时 电价 和 每 个 充 
任务 时 间 必 须 连续 的 假设 下 ， 建 立 电动 汽车 充电 成 本 最 小 模 


ER 


从 现 有 的 研究 文献 来 看 ， 充 电 设施 选 址 问题 都 是 在 经 典 选 
址 模型 上 加 入 里 程 约束 条 件 ， 且 大 多 文献 都 假设 用 户 的 充电 需 
求 是 确定 已 知 的 。 然 而 ， 在 实际 问题 中 ， 先 有 充电 设施 ， 才 会 
有 用 户 来 寻求 充电 服务 ， 而 用 户 需 求 受 到 众多 不 确定 因素 的 影 
响 ， 如 电池 技术 的 发 展 、 交 通路 网 状况 、 用 户 驾 驶 习惯 、 用 户 
心理 因素 、 电 动 汽车 3 此 ， 在 充电 设施 选 址 问题 中 
考虑 需求 的 不 确定 性 至 关 重 要 。 目 前 ， 已 经 有 学 者 将 需求 不 确 
定性 考虑 进 充电 设施 选 址 问题 中 。 Yang 等 人 [3 假设 电动 汽车 用 
户 的 行驶 路 径 已 知 ， 提 出 路 径流 量 不 确定 的 换 电站 选 址 库存 问 
题 。 论文 假设 用 户 流量 是 随机 变量 , 用 概率 约束 作为 约束 条 件 。 
Mak 等 人 04 研 究 基于 充电 需求 不 确定 的 换 电 站 选 址 鲁 棒 优 化 
模型 。Hosseini 等 人 05 研 究 基于 不 确定 性 的 两 阶段 随机 能 源 补 
充 设 施 选 址 问题 。 

现 有 考虑 不 确定 因素 的 充电 设施 选 址 问题 ， 大 多 假设 不 确 
定 参数 概率 分 布 已 知 ， 或 者 不 确定 事件 可 能 发 生 的 情景 已 知 ， 
目标 是 最 优化 期 望 ， 具 有 一 定 的 局 限 性 。 本 文 假设 不 确定 需求 
概率 分 布 未 知 ， 用 盒子 (box) 表示 不 确定 参数 的 取 值 集合 ， 利 
用 基 约 束 鲁 棒 优 化 方法 研究 基于 需求 不 确定 的 换 电站 选 址 库存 
问题 。 首 先 介绍 Bertsimas 和 Simllg 提 出 的 基 约 束 鲁 棒 方 法 ,该 
方法 不 仅 保留 了 确定 问题 的 线性 规划 框架 ， 还 能 够 通过 改变 不 
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确定 预算 来 调节 解 的 鲁 棒 性 。 在 需求 确定 的 换 电站 选 址 模型 基 
出 上 ， 假 设 需求 取 值 于 有 界 对 称 区 间 ， 提 出 基 约 束 鲁 棒 选 址 模 
型 ， 并 将 该 模型 转化 为 等 价 的 线性 规划 。 为 了 分 析 解 的 鲁 棒 水 
平 对 换 电站 选 址 成 本 的 影响 ， 本 文 最 后 应 用 
络 作为 数值 算 例 ， 说 明 解 的 鲁 棒 水 平 对 电动 汽车 换 电站 选 址 成 
本 的 影响 ， 分 析 解 的 鲁 棒 水 平 与 选 址 成 本 之 间 的 权衡 关系 ， 并 
流量 波动 和 续航 里 程 两 个 参数 作 灵 敏 性 分 析 。 


Nguyen-Dupius 网 
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1 不 确定 需求 下 的 换 电站 选 址 模型 


1.1 鲁 棒 优 化 

鲁 棒 优 化 经 Soyster. Ben-Tal, Nemirovski, El-Ghaoui 52?! 
等 的 发 展 ， 被 研究 人 员 广 泛 应 用 于 各 个 领域 254。 本 文采 用 文 
献 [16] 提 出 的 基 约 束 鲁 棒 方法 。 考 虑 如 下 线性 规划 问题 : 


max c'x 


s.t. Ax«b 

l<x<u 
其 中 : xeX ccR" ERRE, pe gn 是 右 端 项 向 量 ， ce R" Æ 
标 函 数 系数 向 量 ，4 < Rw 是 约束 和 矩阵， 其 元 素 是 w 。 对 于 


JERE A HJIT i SJ RRIT i 中 不 确定 参数 的 下 标 集合 ， 


a,c[a,-à,a,*à,]. fa, 表示 不 确定 参数 ，, 表示 均值 或 


.= 
i, y, ef] 
a 


常规 值 ，6, 表示 最 大 绝对 偏离 值 。 令 n= 


且 均 匀 对 称 。 
记 不 确定 预算 为 T, (iel[l,m] )， 


J ={j|6; >0} ， 则 基 约 束 鲁 棒 模型 为 


n 2 à; [xm : 2m; <T, 
S.t. Ya, + max 4 je; je; <b Vi 
ji 0xm,XljeJ, 
l, <x, <u, Vj 
引入 对 偶 变量 wc 18, ， 则 上 述 问题 等 价 于 如 下 线 
性 规划 问题 : 


n 
s.t. Fax, ta, > <b, Vi 
rs 


je; 
oi + B; Zày; , Vi, VjeJ, 
=y; SX; Sy, Vj 
L, <x; Su;, Vj 
a, 20, f, 20, y, 20, Vi, Vj eJ, 
1.2 模型 假设 与 问题 描述 
本 文 模型 作 如 下 假设 : 
a) 模型 只 考虑 纯 电 动 汽车 的 换 电 问 题 , 而 不 考虑 混合 动力 
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汽车 既 可 以 加 油 又 可 以 充电 的 情况 ; 

b) 电动 汽车 用 户 日 常 行驶 的 路 径 已 知 , 即 从 出 发 点 到 目的 
地 的 最 短路 径 。 
c) 电动 汽车 在 起 点 的 电量 为 续航 里 程 尺 的 50%， 在 终点 
的 电量 不 低 于 续航 里 程 R 的 5090051, 该 假设 能 保证 


动 汽 车 往 


Age 
KE 


够 驾驶 电动 汽车 行驶 完整 条 路 径 ， 则 电量 耗 完 
时 ， 必 须 有 换 电站 为 用 户 提供 充电 服务 。 如 图 1 所 示 ， 设 电动 
汽车 的 续航 里 程 g-d00km, ， 用 户 的 行驶 路 径 是 
A->B->C->D， 由 于 续航 里 程 的 限制 ， 用 户 必须 在 B fem 


ce 


池 , 也 就 是 说 只 有 在 B 点 建立 换 电站 , 电动 汽车 才能 继续 行驶 。 
50km B 60km Cy 40km 
图 1 一 条 路 径 了 


对 于 路 径 g ， 首 先 构 造 一 个 扩张 网 络 ， 记 为 G=( 思 ,人 ) ， 


其 中 ys 是 路 径 q TRER, Ñ 
扩张 网 络 G 边 集合 。 d i) 表示 节点 i 和 j LITER 
ord,(i) 表示 节点 i 在 路 径 上 的 顺序 标号 ， 如 图 
通过 以 下 四 个 步骤 构造 路 径 了 


A—B—CoD:» ord, (B)=2. 
的 扩张 网 络 : 
a) 增加 起 始点 9 CAZ 


前 ) 和 


MR. 〈 pp 之 后 )， 并 连 


Bis. 


A 和 4、 D. ER 


1}, Az (s A B,C,D,r) o FEL (s, A) 
M (D, r) 加 入 到 ho 中; 

b) VieN* d (Ai) $5 (si) eA ， 则 路 径 4 ndn 
(s,B); 

c) vie Nd (D) «7 (ir) eh ， 则 路 径 g 增加 边 


(C.t): 
d) 


Vi, j € N° : (ord (i) « ord, Cj)) ^(d, (i, D) < R) — (i, j) e A » 则 路 


ft q 增加 边 (4,B)、(B,C)、(B,D) 和 (C,D) 。 

通过 以 上 四 步 ,路 径 q 的 扩张 网 络 如 图 2 所 示 。 从 起 点 s 到 
终点 1 的 任意 一 条 路 径 均 对 应 一 个 可 行 的 换 电 站 点 组 合 ， 如 路 
有 效 换 电站 点 组 合 为 {B,D} o 


径 y* >B > D >t XMP 


on iy 50km -e 60km 6 40km 
AN Sm AN bkm A o 4- 


图 2 路 径 g 的 扩张 网 络 
本 文 研究 的 问题 是 : 已 知 多 个 用 户 行驶 路 径 的 前 提 下 ， 如 
何 选择 位 置 建立 换 电 站 以 及 确定 每 个 站 点 的 电池 库存 量 ， 以 此 


来 满足 用 户 的 充电 需求 ， 即 通过 
径 。 


在 换 电 站 处 换 电池 来 行驶 完 路 


集合 : 


0 一 路 径 的 集合 ，g e O ; 


参数 : 
c, 一 在 节点 ; 处 建立 换 电站 的 
单位 电池 库存 成 本 ， 
f, P q 上 电动 车 用 户 的 流量 ， 
变量 ， 

[L 在 节点 ;处 建立 换 
> 
7 一 在 节点 ;处 电池 的 库存 量 ， 


| 


h, — Wl 5i i 


By 


ie 


划一 路 径 g 的 扩张 网 络 hy 上 边 (i, 


定 成 本 ， 


ieN: 


先 将 文中 用 到 的 集合 、 参 数 和 变量 表示 如 下 : 


ieN: 


qeQ: 


， ieN: 


N: 


jJ) WWE. Pie 


, 内， 路 径 qg 上 的 车 辆 在 节点 i 处 换 电 
qeQ, ieN’. 
则 
需求 确定 的 换 电 站 选 址 库存 模型 为 
(DM) min $ cy; + A1, (1) 
ieN ieN 
l, i-2s 
st © x- MX w=1-l, iet ,vgqeQ,vieN’ 
(EG je} {(PDs4a9} 0. ist 
Q) 
XxLEm, Vi e N, Vq EQ, (3) 
DUUM 
zt€y,VieN,vVqeQ (4) 
X fu SI, VieN (5) 
4€0; 
x20, V(i, j)e Â’ (6) 
120, VieN (7) 
zí e(061), VgeQ. Vie N* (8) 
y, e {0,1}, Vie N* (9) 


i 


x 
Zq D 车 流 经 过 ys 中 的 节 


式 (4) 表 示 只 有 在 节点 处 建立 换 电站 , 经 过 该 节点 的 车 辆 才能 在 


入 点 处 换 
点 的 车 辆 流量 。 
13 ”和 鲁 棒 选 址 模型 


B. 约束 式 (5) 表 示 换 电站 点 的 
约束 式 (6)~(9) 是 对 决策 变量 的 限制 。 


标 函 数 式 (1) 表 示 最 小 化 总 成 本 之 和 , 包括 换 电站 建站 成 
库存 成 本 。 约 束 式 (2) 表 示 流 量 


均衡 ， 即 每 条 路 径 上 的 
的 地 的 行驶 。 约 束 式 (3) 表 示 路 
点 时 ， 在 该 节点 处 换 电 。 约 束 


电池 库存 不 低 于 


换 电站 的 选 址 问题 中 存在 很 多 不 确定 因素 ,如 用 户 的 需求 ， 
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即 本 问题 中 路 径 上 用 户 的 车 辆 流量 。 而 且 决 策 者 很 难 准确 获得 
需求 的 概率 分 布 信息 ， 本 节 在 上 述 确定 模型 的 基础 上 ， 假 设 路 


径 上 用 户 的 车 辆 流量 是 随机 变量 ， 建 立 换 电站 选 址 问题 的 鲁 棒 
模型 
在 换 电 站 和 鲁 棒 选 址 问题 中 ， 设 路 径流 量 


f. e| f, - ff, AE 其中 表示 路 径流 量 的 均值 或 常规 值 ， 


所 表示 路 径流 量 的 绝对 偏离 值 。 令 |Q|=n(vieN)， 即 路 径 条 


数 为 n 条 。 在 约束 (5) 中 ， 对 vie N ， 引 入 不 确定 预算 T, W 
值 在 区 间 [0,n] 上 。 参 数 的 作用 是 调节 解 的 鲁 棒 水 平 ， 控 制 
解 的 保守 性 。 直观 地 讲 , 对 于 节点 i， 不 是 所 有 的 路 径流 量 f 
会 变动 ， 该 鲁 棒 方法 的 目标 是 当 最 多 有 | 工 | 条 路 径 上 的 流量 在 


E 


其 区 间 内 发 生变 化 ， 一 条 路 径 上 的 流量 波动 为 (r, -LT, f, 


其 他 路 径 上 的 流量 者 

束 条 件 式 (5)。 
需求 不 确定 的 换 电站 和 鲁 棒 选 址 模型 为 

min >》 cy +9 hI, 


ieN ieN 


TN 


取 均 值 时 ， 重 棒 解 仍 是 可 行 解 ， 即 满足 约 


(RM) 


st »AdtRGIOSL WeN (10) 


4€Q, 
式 (2)~(4)(6)~(9) 
其 中 : 


g GL) jrj 


max 
{SU{t}SSQ,, |S|-| T; J. teg; Sj 


(1) 
gc D) Wk 7g W& d? B M. 23 p, REAN, 
gG,T)- sisca Fr) Iia] ! 


命题 1 和 鲁 棒 选 址 问题 (RM) 与 混合 整数 规划 问题 (RM*) 等 


ffr. 
(RM*) min >》 cy +9 A1, 
ieN ieN 
s.t. 式 (2)~(4)(6)~(9) 
> fetar, +94 <I, vieN a2) 
qeQ, qeQ, 
a+, Rf VieN,vqeQ, (13) 
ww>0 VieN (14) 
4,20 VieN,VqeQ, (15) 
证 明 引入 辅助 变量 4 (VsN)， 则 函数 g(r) 可 表示 


为 如 下 线性 规划 问题 : 


(Pi) gG;D)- b> fd: 3 4 ET,0<h, sl,vge o] 
geQ; geQ; 


问题 (P) 的 对 侦 问 题 为 


(DPi) 


s.t. f zi 
a, +A Z fazi NVqeQ, 


a; z0 
4,20 VqeQ, 
问题 (2D) 可 行 且 有 界 ， 由 对 偶 定理 可 知 ， 对 偶 问 题 (DPD 可 
行 且 有 界 ， 则 原 问 题 与 对 偶 问 题 的 目标 函数 值 一致 。 由 此 ， 约 
束 (10) 可 转化 为 以 下 线性 约束 : 


Pfr tolt+ Dh sl VieN 
4€Q, 4€Q; 


a +A Z fazi VieN,vqeQ, 


& z0 VieN 
A,20 VieN,VqeQ, 

和 鲁 棱 选 址 问题 (RM) 与 问题 (RM*) 等 价 , 即 可 得 命题 1 .证 毕 。 

与 确定 模型 (DM) 一 样 ， 鲁 棒 模 型 (RM*) 也 是 线性 规划 。 模 
型 (RM) 可 以 通过 现 有 的 线性 规划 软件 求解 。 由 于 最 优 解 较 容易 
求 得 ， 该 方法 在 实际 问题 中 更 适用 。 

4» pr RI p? AERE RE BUR COM) RO 8 HE BC RM H E 
优 解 ， 显 然 : 


aT;- ,1 20 Vie N,VqeQ, 
4€0; 

即 问 题 (RM3 的 可 行 解 也 是 问题 (DM) 的 可 行 解 ， 所 以 ， 不 
难得 出 Pp* <P” o 

S X 表 示 问 题 (RM*) 的 最 优 解 ， 基 约束 鲁 棒 方法 保证 了 在 
约束 式 (5) 中 最 多 有 了 个 参数 f, 变动 时 ，X "是 可 行 解 。 若 有 多 
Tr,^23 f 变动 ， 此 时 可 证 明 ，X* 以 一 定 的 概率 是 可 行 解 。 
对 于 约束 式 (5)， 和 鲁 棒 解 XY 是 可 行 解 的 概率 可 计算 如 下 : 


"Y Lao J € B(nT,)- 


geQ; 


tl - DH iie AD 


其中 n=|Q|，v=(T,+n)/2， u=v-|v]> 

由 于 参数 能够 控制 解 的 保守 性 ， 参 数 工 也 称 为 保护 水 
平 ， 决 策 者 通过 调节 参数 T, 来 控制 鲁 棒 解 违背 约束 条 件 的 概 
率 ， 从 而 控制 解 的 鲁 棒 性 。 反 之 ， 给 定 解 的 鲁 棒 水 平 ， 通 过 式 
(16) 可 以 计算 参数 T, 的 取 值 。 
事实 上 ， 可 以 利用 二 分 法 来 计算 参数 工 ,的 取 值 。 令 SL* 表 
TBMEIOE. n, 表示 最 大 迭代 次 数 ， 0 表示 误差 容忍 度 。 具 
体 步骤 如 下 所 示 : 


(16) 


令 Eh 


a) l=0 > 


— 05 


st, =1-[a -ctn D Y Z2 kk» Flo) 


Iz] v |21 


... ChinaXive (FRBTII 
RMA a OX, «Wi EGdUDGAX Ub E Cta ARA 
1 100 12 
b) Mia - 7 T, < Mid Bo 
SL, -1-|(1- )C(n|v])* Y C(n1])» keksl: 
' [ HEN 2. | ] 4) 200 —(5)—300— ()—400—(7)— 500—158 
9€ SL «SE » l«-Mid» h «hs € SL »SE» E 
> = G T 
ll,» h «Mid: EW er, Bi: 名 
pm yi F2 (1- 
pan «OH k2n o T eT, Bub WU, $ De 
加 
到 步骤 b)。 和 T 
应 用 基 约 束 鲁 棒 方 法 ， 决 策 者 可 以 通过 选择 恰当 的 参数 T， om 
控制 鲁 棒 解 的 保守 程度 。 更 重要 的 是 ,即使 有 多 于 T, 个 参数 f, x 3. Nguyen-Dupius 网 络 
变动 时 ， 鲁 棒 解 仍然 以 一 定 的 概率 可 行 。 此 外 ， 重 棒 模 型 保留 2.2 数值 结果 与 分 析 
了 原 模型 的 线性 特点 ， 容 易 求解 。 对 于 不 确定 优化 问题 ， 鲁 棒 假设 在 每 个 点 ( 原 有 节点 和 边 上 新 增 节点 ) 换 电站 固定 建立 
解 能 够 满足 约束 条 件 的 概率 称 为 解 的 鲁 棒 水 平 。 成 本 和 电池 库存 成 本 相同 。 路 径流 量 的 最 大 偏离 值 与 路 径流 量 
2 ”数值 算 例 与 分 析 REHE, $ f -02xf,. Hc. AKH [150] 上 随机 生成 13 
2.1 算 例 设计 个 节点 权重 ， 在 鲁 棒 水 平民 不 同 取 值 下 ， 求解 模型 (RM")， 最 优 


本 节 利 用 Nguyen 和 DupiusP231 网 络 图 作为 算 例 ， 如 图 3 所 解 如 表 1 所 示 。 然 后 ， 在 区 间 [1,100] 上 随机 生成 13 个 节点 权 
示 ，Nguyen-Dupius 网 络 图 被 已 被 众多 学 者 用 来 研究 网 络 设施 E, 求解 模型 (RM")， 最 优 解 如 表 2 所 示 。 从 表 1 和 表 2 可 以 看 
选 址 与 交通 问题 。Kuby 和 Lim 指出 当 节 点 之 间 的 距离 较 大 出 当 电 动 汽 车 的 续航 里 程 和 模型 的 鲁 棒 保护 水 平 取 值 不 同时 ， 
时 ， 仅 在 节点 处 建立 换 电站 并 不 能 满足 用 户 的 充电 需求 ， 这 时 利用 CPLEX 求解 模型 (RM) 耗 时 不 超过 30 秒 ,可 见 利 用 CPLEX 
需要 在 网 络 的 边 上 增加 备 选 换 电 站 点 。Kuby 和 Lim j ti — 求解 模型 (RM') 十 分 有 效 。 此 外 ， 当 续航 里 程 相同 时 ,保护 水 平 
种 在 网 络 边 上 增加 备 选 点 的 方法 。 本 文采 用 如 下 方法 : S i, j) FT 越 低 ， 选 址 成 本 越 低 ;保护 水 平 下 越 高 ， 选 址 成 本 越 高 。 


T 


表示 网 络 中 边 (i, py 上 增加 的 节点 , 新 增 的 备 选 点 均匀 分 布 在 边 表 1 节点 权重 在 区 间 [1,50] 上 的 最 优 解 
E: 续航 里 程 R 保护 水 平 工 总 成 本 计算 时 间 ( 秒 ) 
jq 1 19863.6 11.89 
n [eme -1, length(ij)»r 
1(i, j) = r 2 20156.5 14.25 
0, length(i, j) < r " 3 23549. 15.85 
其 中 : length, j) 表示 边 (i, j) 的 长 度 。 该 方法 避免 了 边 长 大 4 25364.1 19.44 
于 7 的 边 上 没有 备 选 换 电 站 节点 的 情况 。 1 12036.8 10.58 
模型 (RM”) 由 AIMMS 3.9 生成 、 CPLEX 12.6 求解 ， 在 2 11548.5 12.55 
100 
CPU 为 Intel Core i7-4710HQ 2.50 GHz, 内 存 为 8.00 GB 的 计算 3 13649.5 19.65 
机 上 进行 实验 。 图 3 表示 13 个 节点 的 网 络 , 每 个 节点 既是 起 始 4 14089.6 23.56 
点 (origins), 也 是 终止 点 (destinations), 一 共有 13x(13-D 二 78 个 ! 8953.6 15.75 
2 2 9015.8 15.64 
LUE : , E 150 
OD 对 ，OD 对 之 间 的 最 短路 即 为 一 条 路 径 。OD 对 之 间 的 流量 3 9987.3 18.62 
(每 条 路 径 的 流量 ) 由 重力 模型 求 得 ， 即 : 4 10034.9 21.77 
f, 三 W;W; x1.5 1 7542.5 16.89 
ij 2 8012.6 19.23 
其 中 : W; ^ Wi; 分 别 表示 节点 i 和 J 的 权重 ， d; 表示 i^ J 之 200 3 8923.4 20.49 
间 的 最 短 距 离 。 4 9232.3 21.67 
1 6534.5 13.54 
2 6926.4 20.49 
250 
3 7136.8 24.09 


4 7582.6 18.91 


".00058v1 


01807 


N 


录用 入 a #, eea a z UA RM 
1 4583.4 18.42 增加 与 解 的 鲁 棒 水 平 之 间 的 关系 。 从 图 5 可 以 看 出 ， 换 电站 选 
2 4956.2 20.54 址 成 本 的 变化 随 着 鲁 棒 水 平 的 提高 变化 越 来 越剧 烈 。 决 策 者 可 
bi 3 5189.4 2309 以 根据 选 址 成 本 的 增加 来 决定 解 的 鲁 棒 水 平 ， 或 者 根据 需要 达 
4 5536.9 24.89 到 的 鲁 棒 水 平 来 计算 成 本 的 增加 量 。 
R2 节点 权重 在 区 间 [1,100] 上 的 最 优 解 表 3 解 的 鲁 棒 性 与 选 址 成 本 变动 
续航 里 程 R 保护 水 平 工 总 成 本 计算 时 间 ( 秒 ) 保护 水 平 。 解 的 鲁 棒 水 平 (%) RE 成 本 变动 (%) 
1 30782.1 15.89 0 45.32 19253.5 0.00 
2 31898.2 16.45 1 49.84 19920.4 0.03 
i 3 326942 18.65 2 54.35 20094.5 0.04 
4 34124.8 19.52 3 58.76 20699.6 0.08 
1 28603.5 12.98 4 63.16 21156.5 0.10 
2 30493.5 1822 5 67.24 22036.9 0.14 
= 3 32841.7 19.65 6 71.33 22986.5 0.19 
4 34510.5 23.56 7 74.92 23936.4 0.24 
1 24369.4 15.75 78.52 24936.2 0.30 
2 25393.5 18.65 9 81.53 25698.1 0.33 
m 3 272637 2047 10 84.54 26936.4 0.40 
4 28065.1 22.56 11 86.93 27956.3 0.45 
1 19863.6 18.92 12 89.33 29068.3 0.51 
2 20156.5 22 58 13 91.14 30934.5 0.61 
zd 3 23549.7 24.92 14 92.95 31865.4 0.66 
4 25364.1 28.47 15 94.25 33053.3 0.72 
1 12036.8 1951 16 95.55 35000.3 0.82 
2 1548.5 2049 17 96.43 36009.5 0.87 
T 3 13649.5 25.19 18 97.31 38323.3 0.99 
4 14089.6 25.81 19 97.88 39523.5 1.05 
i 9534.1 15.86 20 98.46 40690.4 1.11 
2 96974 2461 21 98.80 43909.3 1.28 
m 3 9862.9 28.43 22 99.15 45000.3 1.34 
4 10664.8 29.57 DEN | z 
为 了 说 明 解 的 鲁 棒 水 平 对 换 电 站 选 址 成 本 的 影响 ， 在 区 间 M ; 
[1.100] 上 随机 13 个 车 点 权重 , $ f =0.18x f> n-|g|-78; z 4 E 
a5. : 
BUR A ROO) TESEHUR TI p (CF EEHEJOR. BEKE p 取 不 -E s 
同 值 , 求解 模型 (RM"), 结果 如 表 3 所 示 。 表 3 最 后 一 列表 示 选 "i 
址 成 本 的 变动 ， 也 就 是 保护 水 平 取 不 同 值 时 ， 和 鲁 棒 模型 最 优 成 
本 相对 于 确定 模型 最 优 成 本 的 增加 。 保 护 水 平 与 最 优 总 成 本 之 e DONEC ONE EE 
间 的 关系 如 图 4 所 示 。 由 图 4 可 以 看 出 ， 当 T, =0 时 ， 即 所 有 HKF 
路 径 上 的 用 户 流量 都 取 均值 时 ， 总 成 本 为 19253.5， 此 时 解 的 图 4 解 的 鲁 棒 水 平 与 换 电 站 选 址 总 成 本 
鲁 棒 水 平 为 45.3%， 即 换 电 站 的 电池 库存 能 够 满足 电动 汽车 用 
户 需 求 的 概率 是 45.3%。 随 着 保护 水 平 的 增加 ， 总 成 本 也 增加 。 
当 解 的 鲁 棒 水 平 较 低 时 ,曲线 较 平缓 , 表明 适当 增加 选 址 成 本 ， 
鲁 棒 水 平 就 会 有 较 大 的 提升 。 当 和 鲁 棒 水 平 逐 渐 增 加 时 ， 曲 线 变 


得 陡峭 ， 表 明 通 过 增加 选 址 成 本 的 方式 较 难 提高 鲁 棒 水 平 。 


显然 ， 解 的 鲁 棒 性 是 以 提高 选 址 
量 与 是 否 值得 是 决策 者 要 考虑 的 问题 。 


asp 
m 


成 本 为 代价 的 ， 成 本 的 增 


图 5 表示 选 址 成 本 的 
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ak 
a 


换 电 站 成 本 的 增加 


0.5r 


0.5 


07 08 
解 的 鲁 棒 水 平 


0.6 


图 5 R 


电站 选 址 成 本 增加 与 解 的 鲁 棒 水 平 


23 ”灵敏 性 分 析 


0.9 1 


假设 庆 -o， 对 路 径 用 户 流量 波动 ， 即 o 的 值 ， 作 灵敏 


性 分 析 ， 得 出 在 不 同 鲁 棒 水 平 


， 换 电站 选 址 总 成 本 随 流 量 波 


á x104 
小 一 建站 成 本 
N 电池 库存 成 本 
ab exc --- dbi 
TER 
N 
Wc 
2 * 

1r 

ot ， , i i i J 

250 300 350 400 450 500 550 
电池 续航 流程 

图 7 续航 里 程 的 灵敏 性 分 析 
3 ”结束 语 


电动 汽车 发 展 的 关键 ， 换 电站 由 于 
时 高 电网 负荷 率 、 可 对 电池 进行 慢 


充电 设施 的 合理 布局 是 
有 具 有 前 峰 填 谷 、 电 网 谐 波 、 


动 的 变化 。 当 o= 0.13、0.18 $00.23 时 ， 解 的 鲁 棱 水 平 与 选 址 


总 成 本 之 间 的 变化 关系 如 图 6 所 示 ， 三 条 
相同 : 解 的 鲁 棒 水 平 较 低 时 ， 上 
成 本 ， 可 以 较 大 程度 地 提高 鲁 棒 水 平 。 随 着 鲁 棒 水 平 的 提高 ， 
线 较 陡 峭 , 表明 即使 提高 总 成 本 , 也 很 难 提高 解 的 鲁 棒 水 平 。 


H 


此 外 ， 同一 鲁 棒 水 平 F, H 


越 敏 感 。 


随 着 鲁 棒 水 平 的 提高 ， 三 个 成 本 
随 着 鲁 棒 水 平 的 提高 ， 换 电站 选 


对 电动 汽车 的 绪 


线 的 增长 
线 较 平 缓 ， 表 明 适 当 增 加 选 址 


航 里 程 ， 即 R 的 值 ， 作 灵敏 性 分 析 。 图 7 


趋势 基本 


门 流量 波动 越 大 ， 选 址 总 成 本 越 高 。 
间 的 差距 逐渐 变 大 ， 这 表明 
址 总 成 本 对 路 径流 量 波动 越 来 


表示 电池 的 续航 里 程 对 换 电 站 建站 成 本 、 电 池 库 存 成 本 和 选 址 


总 成 本 的 影响 。 当 续航 是 


8803 降 


EFEM 250 增加 到 500 时 ， 建 站 成 本 从 
到 4601(44.7%), 电池 库存 成 本 从 29632 降 


氏 到 10032 


(66.1%)， 总 成 本 从 38435 降低 到 14633 (61.9%)。 灵 敏 性 分 


电池 库存 成 本 和 选 址 总 成 本 随 续航 里 程 变化 较 大 。 
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图 6 用 户 流量 波动 的 灵敏 性 分 析 


析 的 结果 表明 续航 里 程 的 改变 对 换 电 站 选 址 策略 影响 较 大 ， 且 


[si 


充 延 长 电池 寿命 等 优点 ， 备 受 政府 和 企业 的 推 尝 。 用 户 的 充电 
电站 选 址 模型 的 重要 输入 参数 ， 而 用 户 的 需求 受到 电 
动 汽 车 的 保有 量 、 技 术 的 发 展 、 交 通 网 络 的 状况 等 众多 因素 的 
姑 而 是 不 确定 的 。 本文 假设 不 确定 需求 的 概率 分 布 未 知 ， 
利用 基 约 束 鲁 棒 优化 方法 ， 研 究 基 于 需求 不 确定 换 电 站 选 址 库 
存 鲁 棒 优化 模型 .该 方法 不 仅 保留 了 确定 问题 的 线性 规划 框架 ， 
还 能 够 通过 改变 不 确定 预算 来 调节 解 的 鲁 棒 水 平 。 在 换 
址 确定 模型 基础 上 ， 假 设 需求 取 值 于 有 界 对 称 区 间 上 ， 
约束 鲁 棒 选 址 模型 ， 并 将 该 模型 化 为 等 价 的 线性 规划 。 为 了 分 
析 解 的 鲁 棒 水 平 对 总 选 址 成 本 的 影响 ， 本 文 最 后 应 用 
Nguyen-Dupius 网 络 作为 数值 算 例 ， 算 例 分 析 结 果 说 明 在 不 
的 鲁 棒 水 平 , 换 电 站 选 址 总 成 本 不 同 , 当 解 的 鲁 棒 水 平 较 低 时 ， 
适当 增加 成 本 ， 可 以 较 大 程度 地 提高 解 的 鲁 棒 水 平 。 此 儿 
对 影响 模型 的 重要 参数 作 了 灵敏 性 分 析 。 分 析 结 果 表 明 ， 随 着 
鲁 棒 水 平 的 提高 ， 换 电站 选 址 总 成 本 对 路 径流 量 的 波动 越 来 越 
敏感 ， 续 航 里 程 的 改变 对 换 电站 选 址 策略 影响 较 大 ， 且 电池 库 
存 成 本 和 选 址 总 成 本 随 续航 里 程 变化 较 大 。 本 文 的 研究 以 为 政 
府 部 门 和 电动 汽车 充电 服务 商 决策 电动 汽车 充电 设施 的 位 置 和 
数量 提供 科学 依据 。 
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